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1.  Introduction 


On  atom  chips,  coherent  atomic  waveguides  are  important  components  needed  for  the  construction 
of  certain  types  of  matter  wave  interferometers  that  use  cold  magnetic  atoms.  Results  concerning 
the  use  of  the  time  reversal  symmetry  (Tinkham,  1964)  to  solve  the  quantum  behavior  of  a  spin  one- 
half  particle  guided  in  a  type  of  atom  waveguide  often  used  on  atom  chips  are  presented  here. 

The  Hamiltonian  for  the  atom  guide  is  symmetric  under  time  inversion,  if  the  reversal  of  the 
external  magnetic  fields  is  explicitly  included  in  the  reversal  process.  As  described  in  a  previous 
report  (Golding,  2010a),  the  Frobenius  power  series  method  can  be  used  to  derive  a  pair  of  power 
series  solutions  that  are  well  behaved  at  the  coordinate  origin.  These  two  solutions  are  related  to  one 
another  by  the  time  reversal  symmetry  of  the  system.  Two  other  solutions  that  are  not  finite  at  the 
origin  are  rejected  for  physical  reasons  in  this  waveguide  problem. 

These  series  solutions  are  useful  as  they  allow  the  basic  symmetry  properties  of  the  transverse 
modes  of  the  guided  atoms  to  be  studied  analytically.  In  addition,  these  series  can  be  readily 
summed  for  a  large  range  of  radial  distances  using  arbitrary  precision  arithmetic  techniques  so  that 
the  behavior  of  guided  modes  can  be  studied  numerically. 

In  order  to  form  eigenstates  and  calculate  eigenvalues,  the  pair  of  physically  acceptable  solutions 
derived  in  this  way  must  be  combined  to  satisfy  some  type  of  boundary  condition  relevant  to  the 
problem  of  interest.  Others  have  essentially  used  an  open  boundary  condition  to  describe  a 
quadrupole  guide  of  infinite  transverse  extent.  This  leads  to  an  ambiguity  in  the  eigenstates  that  may 
be  important  in  understanding  the  behavior  of  these  guides  in  practice.  In  our  current  work,  a  hard- 
wall  boundary  condition  is  used,  establishing  a  closed  boundary.  This  forces  both  components  of 
the  spinor  solutions  to  be  zero  at  some  fixed  radius  representing  the  boundary.  When  hard-wall 
boundary  conditions  are  used,  a  complete  system  of  orthogonal  eigenstates  is  found  that  satisfies  the 
Sturm-Liouville  conditions.  Further  work  is  needed  to  understand  the  relationship  between  these 
two  sets  of  boundary  conditions  as  the  limit  of  a  large  radius  for  the  hard  wall  is  taken. 

Since  explicit  series  solutions  are  developed,  an  ordinary  differential  equation  (ODE)  solver  is  not 
needed  to  calculate  the  wavefunctions  or  eigenvalues.  This  is  extremely  useful  as  the  fourth  order 
differential  equations  involved  are  stiff  and  difficult  to  handle  with  normal  ODE  techniques.  Both 
the  energy  spectra  and  the  wavefunctions  can  be  calculated  directly  using  the  series  technique  for 
various  regimes  of  magnetic  field  strength,  allowing  the  direct  study  of  the  assumed  adiabatic 
behavior  (Sukumar,1997)  as  the  longitudinal  bias  field  (Ioffe  field)  is  increased  from  zero  to 
relatively  large  values. 
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2.  The  Atom  Guide  Model 


To  study  the  dynamics  of  an  atom  in  a  magnetic  waveguide,  the  Schrodinger  equation  describing  a 
single  neutral  magnetic  atom  trapped  in  the  magnetic  guiding  field  must  be  solved.  To  minimize  the 
initial  complexities,  the  initial  model  is  chosen  to  be  a  neutral  spin  14  atomic  level.  Two  alkali 
atoms,  40K  and  6Li,  are  fermions  and  have  a  spin  !4  level.  The  important  class  of  bosonic  alkali 
atoms  will  be  treated  by  a  similar  technique  in  a  later  work. 

3.  Hamiltonian  for  a  Magnetic  Atom  in  a  Magnetic  Waveguide 

The  quantum  description  of  an  atom  moving  in  a  spatially  dependent  magnetic  field  must  include 
both  internal  and  external  degrees  of  freedom.  This  means  that  the  momentum,  position,  and  spin 
degrees  of  freedom  must  all  be  treated  as  quantum  operators.  Using  m  for  the  mass  of  the  atom  and 
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M  for  the  magnetic  moment,  the  Hamiltonian  is  written  as  the  sum  of  the  kinetic  energy,  ~  and 

the  interaction  energy  of  the  magnetic  moment  and  the  field  B  .  The  result  is 

2 

H  =  — - M -B(r) 

2m  .  (1) 

The  magnetic  field  B(r  )  is  independent  of  the  z-direction  and  is  taken  to  be  an  ideal  quadrupole 

field  equation  (equation  2)  extending  to  infinity  in  the  x-y  plane  and  uniform  along  z.  An  additional 
uniform  bias  field  B0  is  added  in  the  z-direction  to  help  control  possible  spin-dependent  losses  that 
may  occur  at  the  zero  field  point  at  the  center  of  the  guide.  The  spatial  dependence  of  the  ideal 
quadrupole  field  is  given  by 

B(r)  =  Bl(-xx  +  yy)  +  B0z 

In  the  simplified  model  of  the  atom  described  by  equation  1,  the  magnetic  moment  of  the  atom  is 
just  the  magnetic  moment  of  a  single  outer  electron.  The  atomic  model  used  here  can  be  thought  of 
as  an  alkali  atom  with  a  spin  zero  nucleus  and  a  total  mass  m.  The  magnetic  moment  is  M  =  yS  , 
where  y  is  the  gyromagnetic  ratio  of  the  level  considered  and  S  is  the  spin  angular  momentum  of 
the  atom.  Since  a  spin  14  system  is  being  considered,  the  spin  angular  momentum  is  proportional  to 
-  h 

the  Pauli  matrices,  S  =  —  a  . 

2 

The  Schrodinger  equation  for  the  guided  atom  eigenstate  ^ E  is 
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(3) 


~V2'VA’‘)-rS-B(>!)'rA's)  =  evB(r) 

Zm 

Equation  3  is  made  dimensionless  by  choosing  a  length  scale  such  that  x  — >  he ,  where  the  new  x  is 

h2k2 
2  m 

where  k  =  —  is  the  wavenumber  of  the  optical  field  used  to  cool  the  atom.  This  is  a  rather  arbitrary 

X 

choice  at  this  stage,  but  the  recoil  energy  is  a  convenient  scale  for  atoms  that  have  been  laser 
cooled.  This  is  done  after  dividing  through  by  the  recoil  energy,  pulling  out  the  relevant  factors,  and 

_  kBxX 

defining  the  dimensionless  transverse  field  parameter  U  ~  7  ~  ~  ,  the  longitudinal  field  parameter 

ZbR 

_  hB0  _  E 

bo  ~  7  ,  and  the  scaled  energy  6  -  —  . 

ZhR  Er 


4.  Time  Reversal  Symmetry  of  the  Hamiltonian 


h2  2-  _  frk2 

dimensionless.  Then  —  ^  ^  e  (r)  TZ  ^  (r)  and  the  recoil  energy  is  defined  as  ER 


2m 


2m 


The  Hamiltonian  above  (equation  1)  is  invariant  to  the  operation  of  time  reversal  if  the  external 
field  B  is  reversed  as  well.  This  is  simply  because  when  the  spin,  S  ,  reverses  sign  under  time 
reversal,  B  must  reverse  sign  to  leave  the  Hamiltonian  unchanged.  The  usual  time  reversal 
techniques  do  not  act  on  an  external  magnetic  field.  An  operator  M  is  defined  to  explicitly  reverse 
the  external  magnetic  fields  in  the  calculations,  thus  MBM1  =  —B  .  When  the  dimensionless  field 
parameters  are  present  in  an  expression,  the  field  reversal  operator  simply  changes  the  signs  of  both 
b0  and  bx .  This  operator  is  then  combined  with  the  standard  time  reversal  operator,  T  =  e  ‘S'''T  '  K  , 
where  K  is  the  conjugate  operator  used  in  scalar  time  reversal  problems,  to  form  a  symmetry 
operation  of  the  system  (equation  1).  The  complete  time  reversal  symmetry  of  the  Hamiltonian  is 
then  expressed  as 

(MT)H(MT)1  =H  .  (4) 

Using  this  operator  and  any  given  initial  state  |  w)  an  additional  state  I  ¥ ')  called  the  time  reversed 
state  can  be  fonned  as  | ¥ ')  =  MT \yz)  _ 
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5.  Eigenstates  Common  to  H  and  Az 


As  described  in  Golding  (2010a),  the  equations  for  the  radial  wave  functions  of  the  quadrupole 
guide  are  derived  by  making  use  of  the  angular  symmetry  of  the  quadrupole  magnetic  field.  This 
symmetry  is  expressed  as  the  conservation  of  the  angular  momentum,  A  .sf  -  S_  (Hinds,  2000 

and  2001;  Lesanovsky,  2004;  Bergeman,  1989;  Golding,  2009).  Our  approach  is  to  find 
eigenfunctions  that  are  common  to  both  the  Hamiltonian,  H  and  the  angular  momentum,  A_ .  This 

approach  effectively  separates  the  angular  and  radial  dependence  of  the  problem,  leaving  only  a 
coupled  system  of  radial  equations  (equation  6)  to  solve.  These  equations  are  solved  exactly  using 
both  series  techniques  and  the  time  reversal  symmetry  of  the  Hamiltonian  (equation  1)  to  develop 
the  general  solution  of  the  Schrodinger  equation.  The  general  solution  must  then  be  made  to  satisfy 
appropriate  boundary  conditions  so  that  eigenvalues  and  eigenfunctions  can  be  calculated. 


The  component  fonn  defining  the  radial  spinor  components  for  a  state  with  angular  momentum  pi  is 


R](p)e,,p/1 
R  u  (p)e~ifpl2 


When  the  assumed  fonn  (equation  5)  is  used  in  the  Schrodinger  equation  (equation  3)  the  angular 
dependence  drops  out  leaving  only  the  following  radial  equations  to  be  solved: 


82pR++~dpR+- 

P 


{p  +  \/2f 


R  +(£  +  bll)R_  =  bxpR 


d;R_  +  -dpR- 
p 


R_  +  (s  -  b0  )R_  =  bxpR+ . 


In  the  case  of  a  spin  half  particle  with  integral  orbital  angular  momentum,  the  pair  of  equations  must 

l  3  5 

be  solved  for  the  allowed  half-integral  values  of  the  angular  momentum  ^  ~±2,±^.,±2"'  ■ 


6.  Frobenius  Series  Approach 


When  the  system  of  radial  equations  (equation  6)  is  uncoupled,  a  fourth  order  differential  equation 
results  that  has  a  regular  singularity  at  the  origin.  This  equations  can  be  solved  using  the  Frobenius 
power  series  technique  (Ince,  1956;  Bender,  1978).  The  equation  for  R+  is 


4 


4  d 4  _ 

P  ~r^R+  + 


dp 

( s  A 

2  p'e-  —  +  2p2 
\2 


P 


d2 


dp 


R  + 


2  + 


d 


-2p3b0+(5  +  4p~  +6p^)pj-j-R+  + 

-bfp6  +  [s2  - b2 ) p4  +  f b0  - 2 p2s  +  2 pb0  +  ^s  p2  - y  p2  - 9p - ^  +  //4 


R  =0 


(7) 


A  similar  equation  for  R  can  be  found  simply  by  replacing  R  with  R  and  changing  the  signs  of 
both  ij  and  bo.  This  is  readily  understood  by  inspecting  the  behavior  of  the  coupled  radial  equations 
(equation  6)  under  this  transformation  and  is  due  to  the  time  reversal  symmetry  discussed 
previously. 


The  Frobenius  series  technique  is  applied  to  equation  7  by  assuming  a  modified  power  series 
solution  for  R+  of  the  form 


*.=Xv”' 

-«  ,  (8) 

where  the  coefficients  cn  and  the  index  a  are  unknown  except  that  co  =  1 .  This  is  substituted  into 
equation  7  and  the  resulting  system  is  solved  resulting  in  the  following  recurrence  relations  for  c„, 


cn-6bi2  +  cn- 4  (V  ~  ~  f 2  ( // +  5  / 2  —  w  —  cj)(2s p  -2b0~5s  +  2<577  +  2s cr) 

c  — - - - - - .  (9) 

(2//  +  1  + 2/7 +  2<r)  (2//-7  +  2/7  +  2cr)(2/7  +  5-2/7-2cr)(2/7  +  l-2«-2(T) 

For  n  =  0,  the  denominator  of  equation  9  is  the  indicial  equation  that  is  used  to  determine  the 
allowed  values  of  a  in  equation  8.  The  necessary  conditions  are  just  that  Co  =  1  and  cn  =  0  for 
n  <  0.  The  corresponding  component  R  is  found  by  substituting  equation  8  into  the  first  line  of 
equation  6  and  solving  for  R  producing 

R^Un+a)1  -{n+H2f)p-  .  (10) 

h=o  ^  ’  b^p 

Since  the  allowed  values  of  a  differ  by  integers,  the  maximal  index  is  the  only  index  that  is 
guaranteed  to  produce  a  solution  that  is  finite  at  the  origin  (Ince,  1956).  This  solution  is  called  a 
type  one  solution  and  the  solution  obtained  from  this  by  time  reversal  is  called  a  type  two  solution. 
Type  one  and  type  two  functions  are  denoted  with  a  leading  superscript  such  as  ys  or  R  in  the 
following. 
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In  order  to  use  equation  9  over  the  full  range  of  angular  momenta,  pi,  the  fact  that  the  maximal  index 
cr  is  a  function  of  pi  must  be  considered.  The  maximal  index  a  is  defined  as  a  function  of  pi  by 


<j 


pi>  0 
pi<  0 


(ID 


7.  Time  Reversal  and  the  Two  Solutions 


Equations  8  through  1 1  define  a  single  solution  to  equation  6  for  the  specific  angular  momentum  pi 
that  is  called  the  type  1  solution,  V/( .  The  type  2  solution,  V"“,  is  derived  directly  from  the  type  1 
solution  using  the  time  reversal  symmetry  of  the  system.  In  terms  of  the  radial  wavefunctions,  the 
two  solutions  are  defined  as 


r 

r 


{'R'(b0Ay‘,n) 

'  2K(\,b,)e'”n  ' 


and 


(12) 


The  full  time  reversal  symmetry  MT  acts  on  )//■“  to  give  V“ ,  but  the  fields  and  angular  momenta 
have  to  be  properly  reversed  as  in 

(  In;'/  L  L\yp/2\  f  2n U\J<pl2\ 


MT\^  =eiw 


where 


-'R"_(-\,-by 
K'RZ(-b0,-b,)e-»n  , 


y  (byy 


2R7  (b„,b,)e- 


: <pn 


=  v- 


(13) 


MT  =  -icjMK  =  M 


0  -1 

1  0 


K. 


Thus,  the  time  reversal  operator  creates  a  second  spinor  solution  independent  of  the  first  as  can  be 
demonstrated  by  calculating  the  Wronskian  of  the  solution  pair  in  equation  13.  The  action  of  the 
time  reversal  operator  on  the  spinors  of  different  type  is  summarized  as 


MT  V' (bc,,bl)  =  T  V"(M,)  and 

MT  V (b0,b,)=T  V" =  - V (K A)- 


By  inspection  of  equation  13  the  component  radial  functions  of  the  two  solution  types  are  related  by 
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Application  of  equation  15  to  the  type  one  series  solutions  gives  the  series  solution  for  the  other 
radial  components  in  terms  of  the  single  set  of  coefficients  c„  of  equation  9.  These  series  are 
displayed  in  equation  16: 


XK  (b0 ,bx )  =  p<M)  c„  (< cr ( n ) ,  M,b0 ,bx ) pn 

P 


n= 0 
o-(//)-3  oo 


'R-(boA)  =  iLT — Y,cn(^(p),M^b0,bM(n  +  cr(/j)y -(/u  +  V2f\pn  +  ^^lR^(h0,bl) 


b 


1  n= 0 


bxp 


(16) 


oo 

2R-  {b0A)  =  P^  YaCn(CJ{-p)-P~b 0  ’  ~bl  )  P" 


2K{b0A)- 


P 


n=0 

oo 


b 


1  «=0 


YjCr,(a(-P)’-P~bO-bl) \(n  +  (JirP))  -(A-1/2)“ \pn 


s-b. 


bp 


The  series  expansions  are  used  in  the  general  form  for  the  spinors  (equation  12)  to  produce  a  pair  of 
basic  solutions  to  the  Schrodinger  equation  (equation  3). 


8.  Eigenfunctions  and  Energies 


The  coefficients  cn  (equation  9)  are  functions  of  the  energy  s  .  Some  type  of  boundary  condition 
must  be  used  with  the  radial  functions  (equation  16)  to  detennine  allowed  values  of  the  energy.  One 
possibility  is  to  use  a  hard-wall  boundary  condition  and  force  the  general  spinor  to  be  equal  to  zero 
at  some  radius.  Another  possibility  is  to  assume  there  is  no  hard  boundary  and  to  set  the  general 
solution  equal  to  the  asymptotic  solution  at  some  convenient  distance  from  the  guide  center. 
Experimental  implementations  of  a  quadrupole  guide  will  not  extend  to  infinity  and  some 
understanding  or  model  of  the  physical  boundary  condition  will  be  needed  in  order  to  adequately 
model  the  actual  guide  behavior. 

The  simplest  type  of  boundary  condition  is  the  hard-wall  condition.  This  condition  requires  setting 
the  general  solution  equal  to  zero  at  the  wall.  To  satisfy  this  boundary  condition,  a  linear 
combination  of  type  one  and  type  two  solutions  is  fonned  and  eigenenergies  are  found  such  that 
both  the  general  form  for  both  the  plus  and  minus  components  of  the  spinor  are  simultaneously 
zero. 

The  general  linear  combination  of  the  two  types  of  solution  can  be  fonned  using  a  single  parameter 
0  since  there  will  eventually  be  a  normalization  requirement  on  the  states.  Using  this  parameter,  the 
general  solution  is 
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(17) 


yA  (0)  =  cos(0/2)  y  +  sin  (0/2)  V 

f  (cos(0/2)'<  +sin(0/2)y)e^/2  " 

[(cos(0/2)  lRH  +sin(0/2)  2iA )f?-'*/2 

Assuming  the  boundary  condition  is  independent  of  9 ,  there  are  two  equations  that  must  be  solved 
simultaneously  for  a  common  value  of  0  and  s  .  Using  A,  as  the  radius  of  the  hard  wall,  the 
equations  detennining  the  eigenvalues  are 

cos(6>/2)X',(ff,A)+sin(6>/2)!^'(s,A)  =  0 

(1  o) 

cos(0/2)y  (*,/>o)  +  sin(0/2)2/r  (*,/>„)  =  0 

When  these  equations  (equation  18)  are  solved  a  set  of  orthogonal  spinor  radial  wavefunctions  are 
found  that  may  be  used  for  modeling  the  behavior  of  atom  guides  under  hard-wall  boundary 
conditions.  To  model  open  boundary  conditions  of  infinite  extent,  a  similar  procedure  must  be 
performed  using  the  asymptotic  solutions  of  the  system  (Golding,  2010b). 

When  the  Ioffe  field  is  large  and  the  boundary  conditions  are  open,  discrete  energy  levels  may  exist 
but  the  transverse  modes  are  not  well  defined.  The  lack  of  control  or  consistency  of  the  transverse 
mode  may  lead  to  difficulties  in  the  construction  and  operation  of  beam  splitters  and  other  matter 
wave  components  needed  for  guided  matter  wave  interferometry  and  may  be  one  reason  that  single 
mode  operation  in  these  guides  has  not  yet  been  convincingly  demonstrated.  One  solution  to  this 
problem  would  be  to  develop  a  new  waveguide  design  that  includes  both  a  magnetic  quadrupole 
field  and  some  type  of  conservative  hard-wall  boundary  defined  at  fixed  distance  from  the  guide 
center.  This  may  be  possible  using  a  combination  of  optical  fields  and  miniature  magnetic 
quadrupole  fields,  for  example. 


9.  Conclusions 


The  multiple  solutions  found  in  other  atom  waveguide  work  can  be  easily  derived  using  the  time 
reversal  symmetry  of  a  waveguide  system.  By  organizing  the  series  solutions  using  time  reversal,  a 
more  complete  understanding  of  guide  behavior  can  be  obtained.  For  example,  at  large  Ioffe  fields, 
the  two  solutions  correspond  physically  to  a  bound  mode  and  an  unbound  mode.  The  bound  mode  is 
a  quantum  state  that  is  essentially  attracted  to  the  center  and  its  time  reversed  partner  is  an  unbound 
mode  that  is  repelled  from  the  center.  These  two  modes  must  be  combined  to  meet  prescribed 
boundary  conditions  resulting  in  a  specific  mixture  of  bound  and  unbound  modes  that  can  be 
directly  related  to  the  lifetime  of  waveguide  states. 

This  symmetry  will  be  useful  in  more  complicated  devices  currently  being  developed,  such  as 
magnetic  beam  splitters.  For  example,  if  one  leg  of  a  beam  splitter  favors  the  propagation  of  a 
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bound  state  and  the  other  leg  favors  an  unbound  state,  the  splitter  will  not  work  as  desired.  It  is  also 
hoped  that  the  time  reversal  technique  presented  here  will  be  useful  in  the  analysis  of  atom  guides 
employing  higher  spin  systems.  In  such  systems,  there  will  be  many  possible  modes  involving 
different  spin  states.  Each  state  will  have  a  time  reversed  partner  that  can  easily  be  obtained, 
reducing  the  need  to  develop  many  different  independent  series  solutions.  Time  reversal  techniques 
reduce  the  analytic  and  numeric  work  needed  to  obtain  a  general  solution  as  well  as  enhances  the 
conceptual  understanding  of  atom  guide  systems. 
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